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SUMMARY 

Lagrange’s famous studies of the three-body problem, pre- 
sented in Tisserand’s Traite de Mecanique Celeste , Vol. I, are 
redeveloped in modern mathematical expressions. This new por- 
trayal of a very old problem of celestial mechanics gives rise to 
various considerations which can be very useful for the practical 
computation. In particular, the nine elements of the "reduced" 
three-body problem are defined as a symmetrically constructed 
system of three groups of three elements each, and their differ- 
ential equations are shown. A by-product of this study is a rep- 
resentation of the coefficients of the equation of the fourth degree 
for Lagrange’s quantity p in symmetrical form. 
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ON LAGRANGE’S THEORY OF THE THREE-BODY PROBLEM 


by 

Karl Stumpff* 

Goddard Space Flight Center 


INTRODUCTION 

Lagrange's famous treatment of the three-body problem is still the point of departure 
of many studies on this difficult problem, which has never been solved in its totality. Al- 
though Lagrange's theory is interpretated excellently in its mathematical elegance and 
clarity in Tisserand's textbook of celestial mechanics,** it is intriguing and, for modern 
application, practical to bring forth even more the insight and beauty of this theory by using 
modern mathematical representation. The following considerations are intended mainly to 
serve the end of bringing this classic theory, in new trappings, into the focus of celestial 
mechanical study. Also, in the reconsideration of this old problem, certain previously ig- 
nored questions will be answered. 

Lagrange presented the following postulate: The general three-body problem can be 
solved completely if it is possible to represent as functions of time all those factors which 
are independent of the special selection of the coordinate system. The relative three-body 
problem (i.e., the problem of the motion of two bodies in reference to a third or of all three 
in reference to the common center of gravity) requires twelve integrals, three of which de- 
termine the orientation of the system in space. Thus by Lagrange's theorem only nine in- 
tegrals must be found, for the other three can then be determined easily by a simple quadra- 
ture based upon functions known. 

According to Hesse, we call this problem of seeking the nine "geometric" integrals the 
"reduced" three-body problem.*** It is easy to see that this reduced problem can actually 
be satisfied by nine quantities which are invariable with respect to coordinate transforma- 
tion. The geometric figure formed by the four vectors determining motion the vectors of 
location and velocity of the two bodies with reference to the third — is determined by nine 
quantities which are independent of the coordinate system. These nine quantities are the 
values of these four vectors and the five angles by which the vectors' respective locations 

»NAS-NASA Research Associate; Professor Emeritus, GOttingen University. 

••Tisserand, F., "Traite de Mecanique Celeste ’ Paris: Gauthier-Villars et Fils, 1889, Vol. I. 

•••This is not to be confused with the restricted three-body problem. 
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are fixed. The construction of this figure from these nine elements is not significant in 
itself, since the components can be assembled in several different ways; but of the various 
possibilities there will always be only one that is compatible with the initial conditions 
given. 

These facts, easily understandable geometrically, are reflected in the mathematical 
relationships between the elements of motion. The relative vectors of location of two 
bodies in reference to the third are p , p 2 , and the relative velocity vectors p 1? p 2 ; and 
for convenience let 

Pj “ r i> Vi ~ *V P 2 “ r 3 * P 2 ” r 4* 

The geometrical figure formed by these vectors (discounting the aforementioned additional 
possibilities) is determined by the scalar products possible between them; 

p.b = ( r a r b ) = x . x b + y.y b + z . z b - (i) 

which form the matrix 



Because p ab = p ba , this matrix will be formed from ten different elements among which 
the identity 


Pll 

Pm 

x 1 yi 

z, 0 

P 4 I 

••• P 44 

X 4 y4 

Z 4 0 


holds; therefore only nine of the quantities (Equation 1) which are independent of coordinate 
transformation are arbitrary, whereas the tenth is a function of the others. 


NEW FORM OF THE EQUATIONS OF LAGRANGIAN THEORY 

The beauty of the Lagrange theory can be completely expressed only if we provide 
completely symmetrical formulas. This symmetry cannot be attained if a dominant position 
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is assigned to one of the three bodies, as in the above dis- 
cussion. Therefore we will introduce a reference system in 
which none of the vectors is favored. For this purpose, we 
will call the three masses m 1( m 3 , m 5 (Figure la) or m i f m., 
m k (Figure lb), where the indexes i, j, k traverse the val- 
ues 1, 3, 5 or their cyclic permutations. The relative 
vectors between the masses m. , m k are then labeled 



(a) 


and form the closed train 

r i + r j + r k = °* ( 4 ) 

The same also applies to the relative velocity vectors; if 
we set 



(b) 


l 2» 


or 


r 


a * 


r i 


Figure 1— Relative positions 
of three bodies 


where a = i + l, 0 = j + l, y = k + i traverse the values 2, 4, 6 or their cyclic permuta- 
tions, then it follows from the differentiation of Equation 4 that 


r a + r /3 + r y “ 0 * 

From the six vectors r x , • ■ • r 6 , 36 invariables (r a r b ) - p ab can be formed, but only 21 of 
these are different since p ab = p ba . Twelve relationships exist among these quantities — 
Equation 3 and eleven others — by means of which the five quantities p 15 , p 25 , ■ p 55 and 

the six quantities p 16 , p 66 , ••• p 66 can be reduced to the ten invariables of the matrix 
(Equation 2) using the linear relations of Equations 4 and 5. For instance, 


Pis = ( r i r s) 

( r l> " r l ~ r 3) _P 11 ~ P 13 ' 

(6) 

PS5 = ( r S r 5 ) 

= ( “ r l “ r 3 ) 2 = P U + P 33 + 2P 13 ' 

(7) 


The second of these equations is identical to the cosine theorem of plane trigonometry ap- 
plied to the triangle of the three bodies in Figures la and lb. 
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Equation 3 is not symmetrical since it contains only the invariables with the indexes 
1 to 4; it can be made symmetrical, however, if we use Lagrange’s quantity p. 

Lagrange noted that the remarkable relationship 

2 P = Pm" P23 = P25 “ Pi« = Pj« " P*s ( 8 ) 


exists among the invariables, for 

Pm " P23 = “ (p 1 6 + P12) + (p 2S + P21) = Pis - Pis' 

and 

Pm " P23 = " (P34 + P54) + (P43 + P# 3 ) = P 36 " P4s- 
On the other hand, if we set 

2s j = P 38 + P4 S ' 

2S 3 = P2S + Pl6- 

2s s = Pm + P23' 

from Equations 7 and 8 we obtain 



Pm " s - 

; + 

P 2S = s 3 + p, 

P 36 

: s, + P 


P 23 = s s 


Pl6 = S 3 " P' 

P 4 S = 

Z S 1 ~ P- 

(10) 

If Equation 10 is substituted into Equation 3 for p 14 

and Pj3 . 




Pu 

P 12 P 13 

s s + P 



M = 

P 21 

P 2 2 s 5 “ P 

P 24 

- 0 , 

(ID 


P 31 

S 5 ” P P33 

P34 




S s + p 

P42 P43 

P44 



which is a fourth degree equation for p in which the cubic term is missing: 



p‘ * + 

■ Ap 1 + Bp + C = 

: 0 . 


(12) 
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The coefficients of this equation are not constructed symmetrically however, if they are 
determined by the solution of the determinant (Equation 11). But since p is independent of 
the mass favored (Equation 8), it must be possible to give to the quantities A, B, and c a 
symmetric form based upon the three masses. 

So far no effort has been made to show this in extenso . In order to do so, nine quan- 
tities which form a symmetrical system will be selected from the 21 invariables 
Pn- Pi 2 - " Pss- such as p,,, p 33 , p 55 ; p 22 . p 44 , p 66 ; p 12 , p 34 , p S6 , or the squares of 
the relative distances and velocities, and the scalar products of the three vectors of 
location with their respective velocity vectors. These nine quantities shall be called the 
fundamental invariables of the reduced three-body problem. If these quantities are known 
as functions of time, all others can be derived from them, some with linear relationships 
like Equations 6 and 7, some with the help of the Lagrange quantity p, which, because of 
Equation 12, will be a function of the fundamental invariables. Twelve of the 21 invari- 
ables can therefore be expressed with p and the nine quantities. In actuality, we find 

2 P 13 = 

2p 24 = 

2p 3 s = 

2P 46 = 

2p S i = 

2P 62 = 

Further, 

2s i = 

2S 3 = 

2s 5 = 

and the remaining six invariables are obtained from Equation 10. 

With the help of these relationships, it is also possible to express the coefficients of 
the biquadratic equation (Equation 12) with the nine fundamental quantities in symmetrical 
form. We shall present only the results of this somewhat complicated computation. If, 


P 55 “ P 1 1 ~ P 33 * 


P 66 P 22 P 44 * 


P 1 I P 3 3 ~ P 5 5 


P 22 ~ P 44 ~ P 6 6 


P 3 3 “ P 5 5 “ P 11 


P 44 ” P 66 “ P 2 2 


(13) 




J \2 ^34 ^56 ’ 


P 34 ~ P 5 6 ” P 1 2 


P 5 6 P 12 ” P 3 4 » 


(14) 


J 
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for purposes of brevity, we introduce 

*11 = P11P22 " Pi 2 2 - 

*13 = P11P44 “ P34P12 • 

*15 ~ Pll P 66 ~ P 56 P 12 ’ 

*31 ” P33P22 ~ P 12 P 34 ’ 

*3 3 ~ P33P44 ~ P 3 2 4 ’ 

*35 = P33P66 ‘ P 56 P34 > 

*S1 P 5sP22 P 1 2 P S 6 ’ 

*53 ~ P55P44 ” P34PS6 * 

*5 5 “ P 5S P 66 ~ P 5 2 6 ’ 

X 13 = *13 + *31 ’ 

X 3 5 = *35 + *53 ’ 

X 51 = *51 + *15 ■ 

we obtain 

2A = *„ + ^33 + *55 - ( X 13 + X 3 5 + x 5l) ■ 

P 11 P 33 P 5 5 

P 2 2 P 44 P 6 6 

P 1 2 P 3 4 P 5 6 
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The unhandy necessity of solving Equation 12 can be avoided in all practical applica- 
tions by using the always available integral of the constant impulse moment. Actually, the 
two known integrals of the reduced three-body problem can be expressed through the fund- 
amental invariables andp. We have the energy integral 


h 


1 ft* 2 2 + P 44 + P 6g\ / 1 + 1 1 \ 

2 ' 1 + + ”*/ V-n.yPTT +m a/P73 + * 5 TO/ 


constant ; 


(17) 


and, if 


'Pil - PiiPaa ~ Pi a = 1, 3, 5; a = 2, 4, 6) 

are the diagonal terms of the matrix t p lk (Equation 15) and if 

2<*> jk = A + X jk + 2 p 1 , 


the impulse moment is 


r 



■ in^ i uij 

g 2 = G + m t m 3 m s ( A + 2 P 2 ) = constant. (18) 

In Equation 18 g is the length of the impulse moment vector which stands vertically upon 
the invariable plane, and 

G = hi + h± + ^ + hl_ + Al + _^h. 

m i 2 m 3 m s m 3 m S m 5 m i m i m 3 

If the unit of mass is chosen so that m 1 +m 3 + m 5 = i, it follows from Equation 18 that 

2 p 2 = m 1 m 3 m 5 (g 2 - G) - A (19) 

or, if p 2 and p 4 are eliminated with the help of Equation 19 from Equation 12, 

A 2 1 

= ~ C - *4 m 2 m 2 m 2 ( g 2 - G ) 2 • 


B p 


( 20 ) 
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Since g 2 is known from the initial conditions of each special problem, p can be obtained 
from Equation 19 or, if B is not equal to zero, from Equation 20, 

The differential equations of the relative three-body problem can be written in vecto- 
rial form: 


Mi 


r, + m. R 


(i = 1,3,5; I m. = l) (21) 


where 


and 





y ( 22 ) 


R = r 1 p 1 + r 3 /i 3 + r 5 M 5 


These differential equations are not independent of one another because of the rela- 
tionship saown in Equation 4, Thus it is sufficient to solve two — perhaps those for rj , r 3 
(two vectorial differential equations of the second order which form a twelfth order sys- 
tem). But this system is not symmetrical since one of the three masses, perhaps m 5 , is 
given a favored position. Instead of this, however, the Lagrange theory makes it possible 
to relate the problem’s solution to the integration of nine first order differential equations 
for the nine fundamental invariables, i.e., a completely symmetrical system with no un- 
necessary components. In actuality, if we differentiate the quantities 


Pii = ( r i r i)> P,a = ( r l r a)’ Paa = ( r a r a) !i “ 1. 3, 5; a - i +1) (23) 

with reference to time and eliminate r a = r, by means of Equation 21, then 


1 . 

2 p ii = P ia ■ 


i. - p “p -Lf. ^m.p. 
* i a F aa r; l 


p - ** p . a. + m . p 

2 K aa p iar i i ^ a 


> (24) 


in which we abbreviate 


Pi = ( r i R ) - Pil^l + Pi3^3 + PiS^S 

. = (*. B ) = 


(i = 


Pc 


Paj^l + Pa3^3 + Pa 5 ^5 ' 


a 


1 , 3 , 5 ; 

2 , 4 , 6 ) • 


(25) 
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The formation of the quantities p. and p a from the fundamental invariables becomes even 
clearer if we write 


where 


Pi - Pijfo - Mi) + Pik(^k " Mi) . 

Pa = s k(Mj - Mi) + s i (M k - Mi) + M (M k “ Mj) . 


(26) 


2p : 


r Pk 


2s, 


“ Pi 


ky 


Pi, 






If Equations 24 are integrated (by using, for example, the numerical methods for the given 
initial conditions), the problem is solved in principle and the integrals (Equations 17 and 
18) are available to check the results. The Lagrange's differential equation for p, 

2 P ~ m ,P 35 K " ^s) + m 3P 5 lK " ^l) + m 5 P 13(^l " 


= L " h.) ' 


(27) 


can be added to the system (Equation 24) as a tenth equation, and is remarkable in its sim- 
plicity and symmetry, thereby making possible another thorough check. 

Under certain conditions — especially if one of the three bodies has a considerably 
larger mass than the other two and their orbits can be considered to be disturbed Kepler 
movements in conical sections — it is useful to introduce, in place of the fundamental in- 
variables, derived invariables which are constructed to correspond with those used in the 
theory of the two-body problem.* If we set 




l 



Pii 


3 / 3 


cr. 



p. - 2p i ~ oij , 

e i = “i " ^ ■ 

<5 . = &j. - a. 2 , 

1 1 1 


y (28) 




Stumpff, K., "Calculation of Ephemerides from Initial Values,” NASA Technical Note D-1415, in publication 1962. 
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then these quantities, as derived from Equation 24, satisfy the differential equations 


Mi 


3/^(7. (or r, - r.cr.), 


a. - e. -2 a? + m. — — 
1 1 1 1 Pii 


2o- i (/x i + %) + 2 mi — 


P i 


2p . cr. - 2m . 


i P: : ’ 


r. (a. + 2e) + 2m. ~ f 

i \“i i i p • - * 


= -475.(7, + 2— (p a - PjCT.) 


(i = 1, 3, 5; 

a = i + 1 ) 


r ( 29 ) 


Since only three of the six quantities (Equation 28) are independent, only three of the equa- 
tion systems (Equation 29) need to be considered, i.e., e.. If we also set 


Pi 


e . = e . + m 


i “‘l p. , 


(30) 


we obtain the system: 


where 


Mi = - 3^i . (or = r.a.) 

= e* - 2a -2 ; 

r i (Mi + 2«*) + y, , 


• * _ 

e, - -a 


(31) 


ill . 

r~ (2P a + Pi) 


7i 


(32) 



11 


The y are three functions of the invariables which, if m . is one of the small masses, are 
called perturbations. To show this, we form the equations 


? i = f + r i^i = r i ( £ t - a i) ■ 
-f 'i = f i( e t _cr i J ) + r i( £ t - ^i) 


- r.o-j ^3 (e * - a*) + M .] + 


by differentiating the identity r* = x . 2 + y . 2 + z. 2 . After elimination of o-. and e* - a* from 
these equations, there are three third order differential equations: 


r.r. + 3r . r . + — r 

ii ii ^2 


r . 

i 

T: 2 


r ? y. 

i ' i 


"i ( 2 Pa + Pi) • 


(33) 


whose nine integrals also solve the reduced problem. Because lp a = 1 p. = o, it follows 
from Equation 33 that 



0 , 


an equation whose integral 



constant 


yields the energy law. 


CONCLUSION 

In the case of general three-body motion where the ratios of mass can have arbitrary 
values, we always have to rely on numerical integration methods. We have a choice of in- 
tegrating either the nine first order differential equations (Equations 24 or 31) or the three 
third order differential equations (Equations 33). The preferable method is determined by 
experience and available facilities. The system (Equation 33) would probably be preferred 
for program- conti'olled electronic computers, since it requires only three, rather than 
nine, tables of differences. 
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For instance if the problem is to calculate disturbed planetary orbits, when mi denotes 
the mass of the sun, and m 3 , m s two small planet masses (such as Jupiter, m 3 < lcr 3 and 
Saturn, 3m s < 1CT 4 ), and if the perturbations of the second order can be ignored, we use the 
six equations (Equations 24 or 31) or the two equations (Equation 33) with i “ 3, 5. The 
disturbing functions r 3 , y s , which are compounded with small factors, are then calculated 
with the help of the system’s undisturbed movement which is known from the initial condi- 
tions. This process cannot be used in equations where i = l; the motion of the masses 
m 3 and m 5 relative to one another cannot (as can that of m 3 or m 5 relative to ) be consid- 
ered in first approximation to be a Kepler ellipse. It is unnecessary, on the other hand, to 
solve the equation where i - l, since its solution is implicitly a co-product of the integra- 
tion process leading to the solution of the other equations. 

It should also be mentioned that 


and 


-^( r i 2 ^i) - -2m.p a 





2m iPii (Pa " Pi 0 '!) 


> (34) 


can be derived from Equation 29. These equations are integrable for m. = o and yield the 
known integrals of the two-body theory: 


r 2 p. = -L = constant (energy theorem), 


r 4 $ - p. = constant (area theorem). 

ii i 

In this method of writing the invariables, the formulas (Equation 34) represent four of the 
differential equations of the elemental perturbations used in the theory of the special 
perturbations. 

The task of carrying out the integration of the equations of perturbations according to 
a method which has been earlier described for the solution of the undisturbed problem will 
be treated in a later report. An attempt of this type has already been made for Hill’s lunar 
problem,* a particularly simple variation of the restricted three-body problem. The ex- 
pansion of this result to include the general problem of disturbed orbits will be of consid- 
erably more interest. 


•Stumpff, K., "Remarks on Hill’s Lunar Theory, I and H," NASA Technical Notes D-1450, D-I451, in publication 1962. 
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